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ABSTRACT 

We present high spatial resolution (FWHM ps 0.3-0.5") H (1.6/xm) and K' (2.1/im) 
images of a complete sample (12) of ultraluminous infrared galaxies (ULIGs) chosen 
to have "warm" mid-infrared colors (/25//6O > 0-2) known to be characteristic of 
active galaxies. The extended underlying galaxy is detected in each system at H 
and K 7 , as are tidal features and many of the star forming knots seen at optical 
wavelengths. While some of these knots have a considerable near-infrared excess, we 
show that they are likely to have bolometric luminosities more similar to that of the 
extended starbursts with similar optical morphology seen in less far-infrared luminous 
interacting systems. We find that each ULIG has increasing contributions at long 
wavelengths by a very compact source which we identify as an active galactic nucleus 
(AGN). We show that the optical/near- infrared colors of these putative nuclei are 
much more extreme than the most active starburst IRAS galaxies, yet are identical to 
"far-infrared loud" quasars which are in turn similar to optical quasars combined with 
large quantities of hot (~800 K) dust. Half of the ULIGs have nuclei with dereddened 
near-infrared luminosities comparable to those of QSOs, while the other half have 
dereddened luminosities more similar to Seyferts, although this may be an effect of 
patchy extinction and scattering. 



Subject headings: infrared: galaxies — galaxies: star clusters — galaxies: galaxies: 
active — galaxies: starburst 

Author's note — due to space considerations in this preprint, the figure quality has been 
downgraded. A higher quality preprint can be obtained from http://spider.ipac.caltech.edu/ 
jason or by writing the author directly. 
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1. Introduction 

One of the most important results from the Infrared Astronomical Satellite 1 (IRAS) all-sky 
survey was the discovery of a significant population of galaxies that emit the bulk of their 
luminosity in the far-infrared (e.g. Soifer et al. 1984). Studies of the properties of these "infrared 
galaxies" showed systematic trends coupled to the total far-infrared luminosity; more luminous 
systems were more likely to appear to be merger remnants or interacting pairs, and were more 
likely to possess AGN-like emission line features. A more complete review of the properties of 
luminous infrared galaxies can be found in Sanders & Mirabel (1996). Much recent attention has 
been focused on ultraluminous infrared galaxies (ULIGs), objects with infrared luminosities, L; r 2 , 
greater than 10 12 Lq, which corresponds to the bolometric luminosity of QSOs 3 ( assuming the 
blue luminosity criterion Mb < —22.1, adjusting for our adopted cosmology : Schmidt & Green 
1983). Multiwavelength observations of a complete sample of 10 ULIGs led Sanders et al. (1988a) 
to suggest that these objects might plausibly represent the initial dust- enshrouded stage in the 
evolution of optically selected QSOs, and that the majority, if not all QSOs may begin their lives 
in such an intense infrared phase. 

An important subset of ULIGs are those objects with "warm" mid-infrared colors 
(/25//6O > 0.2). 4 These warm objects, which represent ~ 20-25% of the total population of 
ULIGs discovered by IRAS, appear to represent a critical transition stage in the evolution of 
the larger population of "cool" ULIGs into optical QSOs. Studies of several small but complete 
samples of warm ULIGs have shown that many of these objects have a point-like optical appearance 
on the Palomar Sky Survey and that they exhibit broad (i.e. Seyfert 1) optical emission lines, 
characteristics that have led them to be referred to as "infrared QSOs" (e.g. Low et al. 1988, 
Sanders et al. 1988b). This is not surprising since this criterion has been used previously to select 
AGN from the IRAS Point Source Catalog (de Grijp et al. 1985, 1987). Particularly useful for 
study has been the complete flux-limited sample of 12 warm ULIGs from the survey of Sanders 
et al. (1988b). As the nearest and brightest warm ULIGs, these objects are the most amenable 
to studies at other wavelengths. Other selection techniques (e.g. 12/um flux selection; Spignolio & 
Malkan 1989) have produced similar, albeit smaller samples of ULIGs with AGN spectra. 



lr The Infrared Astronomical Satellite was developed and operated by the US National Aeronautics and Space 
Administration (NASA), the Netherlands Agency for Aerospace Programs (NIVR), and the UK Science and 
Engineering Research Council (SERC). 

2 Li r = L(8-1000/im) is computed using the flux in all four IRAS bands according to the prescription given in 
Perault (1987); see also Sanders & Mirabel (1996). Throughout this paper we use H a — 75 km s _1 Mpc _1 , q a = 0.5 
(unless otherwise noted). 

3 Based on the bolometric conversion i bo i = 16.5 xvh v {S) of Sanders et al.(1989) for PG QSOs. Elvis et al.(1994) 
indicates a value of 11.8 for UVSX QSOs, increasing Mb to -22.5. 

4 The quantities /12, /2s, /60, and /100 represent the IRAS flux densities in Jy at 12 /jm, 25 /Km, 60 /u,m, and 100 /im 
respectively. 
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Surace et al. (1998; Paper I) examined the warm ULIG sample of Sanders (1988b) using 
ffi>T/WFPC2 and found that nearly all were advanced merger remnants with extremely complex 
optical morphologies in their central few kiloparsecs. These central regions were dominated by 
knots of powerful star formation plus one or two compact emission sources whose optical colors, 
luminosities, and physical sizes appeared similar to reddened quasars, and were therefore presumed 
to be putative active galactic nuclei (AGN). Unfortunately, the unknown effects of extinction on 
small spatial scales resulted in considerable uncertainty in the properties of the star-forming knots 
and putative AGN, and optical colors alone could not effectively differentiate between the two 
phenomena. 

Extending wavelength coverage to the near-infrared offers several potential improvements 
over previous studies at optical wavelengths. The most important of these is lowered sensitivity 
to extinction effects; the extinction at 2.1/mi is only one-tenth that at 0.5 fj,m. In these objects, 
which are expected to have some regions of significant dust obscuration, the near-infrared offers a 
much greater opportunity to detect additional emission regions which may be so heavily obscured 
as to be invisible at optical wavelengths. While some recent studies of ULIGs have found that 
their central regions are opaque even at near- infrared wavelengths (e.g. Goldader et al. 1995), 
these studies have concentrated primarily on "cooler" systems which are likely to be more 
heavily obscured. The presence of optical Seyfert lines in the warm ULIG sample would seem to 
indicate that extinction is much less severe than in the cooler ULIGs. Also, determination of the 
near-infrared spectral energy distribution (SED) of the optically observed small-scale structure 
allows a disentanglement of the effects of extinction, star-formation, and AGN activity. Color 
information at longer wavelengths is particularly important since this is where normal stellar 
populations and AGN have very different colors and because thermal emission from hot dust first 
becomes detectable longward of H-band. 

High resolution observations greatly increase point source sensitivity; in the background 
limited case doubling the spatial resolution is twice as effective as doubling the telescope aperture. 
The increased point source sensitivity in turn enhances the detectability of features such as 
star-forming knots. The distributed spatial morphology of these star-forming knots allows 
them to be distinguished from AGN. By identifying these star-forming regions (some of which 
may be hidden from optical detection by dust extinction) and characterizing their SEDs in the 
near-infrared, a more accurate assessment of the contribution of star-formation to the total ULIG 
energy budget may be obtained. Unfortunately, the characteristic size-scale of the optical structure 
is too small (ssO.2-0.4") to be resolved by conventional means even at the best of ground-based 
sites. However, recent advances in adaptive optics (AO) techniques and deconvolution, which 
compensate for local and atmospheric image distortions, allow this kind of resolution to be 
achieved from the ground. We report here the results of such a high spatial resolution near-infrared 
imaging study using a low-order tip/tilt system. 
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2. 



Data 



The data were taken between October, 1995 and May, 1997 with the UH 2.2m telescope on 
Mauna Kea. The observations were made at H (1.6/xm) and K' (2.1/xm). The H filter was chosen 
since it is the longest wavelength filter which is still relatively unaffected by thermal dust emission; 
dust hot enough to emit significantly at this wavelength would be above the dust sublimation 
temperature. The choice of the University of Hawaii K' filter, which is bluer than both the 
K-filter and the 2MASS Kg, was motivated by the lower thermal sky background in the K'-band 
(Wainscoat & Cowie 1992). This improves the detectability of faint features as star-forming knots. 
Throughout this paper we exclusively refer to K' . Comparison to work by other authors is made 
using the conversions (Wainscoat & Cowie 1992): 



A fast tip/tilt f/31 secondary, capable of correcting for image motion at a rate of 10- 100 Hz, 
was used to stabilize the images against both correlated seeing effects such as telescope wind shake 
as well as some atmospheric seeing (Jim 1995). The detector was the QUIRC 1024x1024 camera 
utilizing a HAWAII HgCdTe array (Hodapp et al. 1996). At the f/31 focus, it's 0.06"pixel" 1 scale 
ideally samples the focal plane under diffraction-limited conditions (0.2l"at 2.1//m). Because of 
the camera's relatively large field of view (^61"), it was possible to dither the observations on 
size-scales larger than the detectable galaxies themselves and subtract successive exposures in 
order to correct for the sky background, thus doubling telescope efficiency. The actual pattern 
used depended on the structure of the galaxy. Pixel-pixel response variations were corrected 
by dividing the sky-subtracted images by high S/N normalized dome flats. Based on poisson 
statistics the dome flats are expected to be accurate to greater than 0.1% and the stability of the 
array generally allows flattening with residuals far below 1%; in the near-infrared one is almost 
always more limited by the bright background sky. The dithered frames were registered using 
the centroids of the background stars and galaxies in the frames or in some cases with cross- 
correlation using the extended features in the ULIGs themselves. In both cases the registration 
should be accurate to within « 0.25 pixels. The images were scaled by their exposure times and 
any constant residual background was subtracted. They were then combined using a clipping 
algorithm in IRAF/IMCOMBINE that rejects high-sigma outliers and weights the images by 
their exposure times. Typical total exposure times were 30-120 minutes per target per filter. In 
general, any given target was observed using the same exposure times. Different exposure times 
were necessitated by the brightness of the objects and the seeing conditions at the time. The 
QUIRC camera has a very linear response over nearly its full range; all exposures were well within 
the linear regime. 



K-K' = 0.18 (H - K) 



(1) 



(H - K') = 0.82 (H - K) 



(2) 
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The data were flux calibrated by observing several photometric standards (Elias et al. 1982, 
Casali & Hawarden 1992) throughout the night at the same airmass as the targets. In most cases 
the nights were photometric with a characteristic residual scatter of < 0.05 magnitudes. For 
those where they were not, we calibrated the data using either earlier photometric observations 
at the Palomar 5m (taken with a 58x62 InSb array between October 1988 and March 1992) or 
the aperture photometry of Neugebauer et al. (1987) and Sanders (1988a). The K' filter was 
calibrated using the same infrared standards and the above K—K' correction. 

The point spread function (PSF) was calibrated in most cases by using actual stars in the final 
combined image. Point sources were identified by eye and the DAOPHOT package in IRAF was 
then used to construct a model PSF for each image. This models consists of a theoretical elliptical 
gaussian core and a "correction image" representing the observed variations from the theoretical 
model. In a few cases there were no suitable stars in the vicinity of the target. In these 
PSF was obtained by interleaving observations of the ULIG and nearby bright stars by using the 
same guide star (and by extension, nearly the same degree of image stabilization) for both sets 
of observations. Although the derived PSF is not necessarily identical to that of the ULIG data, 
there is evidence that the atmospheric variability on Mauna Kea is stable on timescales of tens 
of minutes (Ramsay et al. 1992), comparable to the interleave times used. Typical of most AO 
systems, the observed PSF consists of two components: a diffraction-limited core and an extended 
halo approximately 0.5-1" in diameter, equivalent to the uncorrected PSF width (Figure 1). This 
extended halo is due to scattered light in second- order and higher distortion modes, which cannot 
be compensated for by a tip/tilt- type system (Beckers 1993). Although on a small telescope 
(such as the UH 2.2m) tip/tilt is capable of delivering nearly all the gain in resolving power that 
can be achieved with higher-order systems, this uncompensated extended halo can still contain a 
significant fraction (« 30%) of the total flux, thus complicating photometry. All photometry was 
performed using aperture photometry in circular or polygonal apertures, depending on the shape 
of the observed structure. Aperture corrections were derived from the observed PSF. 

In several cases, deconvolution was applied. This was particularly necessary in those cases 
where structure was seen on very small spatial scales (< 0.5" ). Deconvolution is able to increase 
the contrast and detectability of small-scale structure and improve the photometry by eliminating 
confusion resulting from the extended scattered halo. We chose to employ the Richardson-Lucy 
algorithm (R-L) implemented in IRAF. The algorithm was allowed to iterate 25 times. As a 
check on the validity of the deconvolved structure, we compared the deconvolved images with the 
HST optical images, which are of much higher resolution. In general they have similar structure, 
and will be discussed more in §3.1. In all cases photometry was performed on both the raw and 
deconvolved data. 

In the case of one galaxy, IRAS 07598+6508, it has not been possible to obtain complete, high 
quality data at both H & K 7 . Therefore, while images and radial profiles at J and K-band are 
presented, it will be excluded from further discussion. It has a Seyfert 1 optical nucleus, leading 
some authors to describe it as an "infrared quasar" (Low et al. 1988), as well as considerable tidal 
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Radius (arcseconds) 

Fig. 1. — Typical radial profile of the point spread function at K 7 achievable with fast tip/tilt 
guiding. A FWHM of A/D at K'has been denoted by the scale-bar . Note that the central core is 
nearly diffraction-limited in width, but that a considerable amount of power remains in the wings 
due to higher-order scattering modes. 
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debris in the surrounding host galaxy (Boyce et al. 1996). Since it is widely regarded as a QSO, 
it will therefore be considered similar to the three "infrared-loud" QSOs newly observed at H & 
K' and discussed in §3.3.4. 



3. Analysis 

3.1. Near-Infrared Morphology 

ffi> T/WFPC2 imaging indicates that the optical emission of these warm ULIGs is characterized 
by a dominant nucleus which is the most compact (generally unresolved by HST) and luminous 
source at B and which has been identified as a putative AGN, as well as by compact star-forming 
regions (with the exception of 3c273) found predominantly in the inner few kpc (Paper I). In 
most cases these star-forming regions seem to have ages based on their optical colors of less than 
w 5x 10 8 years. The near-infrared, however, traces primarily the late-type stellar population at 
H-band with increasing sensitivity to emission from hot dust at K'-band (Aaronson 1977). Figure 
2 shows the large scale structure observed at H-band for the complete sample of warm ULIGs. 
The near-infrared morphology generally appears to be much simpler than that seen optically. In 
every object, the extended emission of the surrounding galaxy is detected and in several cases 
we can even trace the tidal structure (e.g. Mrk 1014, IRAS 05189-2524, IRAS 08572+3915, Mrk 
231, Pks 1345+12, and Mrk 463). Three of the objects (IRAS 08572+3915, Pks 1345+12, and 
Mrk 463) show clear double galaxy nuclei, a result known from previous optical and near-infrared 
imaging (Sanders 1988a; Carico et al. 1990). Additionally, we can detect the spiral arms and 
star forming knots of the host galaxy of I Zw 1, the star- forming knots in the outer regions of 
Mrk 463 and IRAS 15206+3342, and the extended host of IRAS 12071-0444. However, only in 
a few cases (IRAS 15206+3342, Mrk 231, and Mrk 463w) is there any indication of complex 
near-infrared structure in the immediate vicinity of the nuclei, and even then this structure is 
comparatively faint compared to the nuclei themselves (Figure 3). This smoothness is in part due 
to the lower resolution of the near-infrared images. However, the typical size scale of the optical 
knots is 150 pc, with a knot separation of 300-500 pc. In this case we would expect many of the 
knots to be separated. Typical B—K* colors for a starburst lie in the range 0-4, with B— K 1 ps 
2.5 for a starburst 10 8 -10 9 years in age. Given that the starburst knots detected in Paper I were 
typically me=22-24, then most of the knots should be mx' ~19.5-21.5. The detection limit for 
point sources is typically tok' =22-23; most of the optical knots should therefore be detectable. 
Unfortunately, matters are complicated by the more limited spatial resolution in the near-infrared, 
which limits the dynamic range in the vicinity of very strong point sources, and also by confusion 
due to the much stronger emission in the near-infrared from the underlying host galaxy. The 
former may not actually be very significant, given the nearly unresolved nature of the K-band 
peaks (see below). 

Despite the more limited spatial resolution in the near-infrared, the star-forming knots seen 
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Fig. 2. — The large-scale structure observed at H. The data have been adaptively smoothed, a 
process which applies greater smoothing to regions with low signal- to- noise in order to increase 
contrast. The data are displayed such that the gray- scale part of the image illustrates low surface 
brightness features, while the white logarithmic contours show the high surface brightness structure. 
Tick marks are 0.5 "apart, with major ticks every 4 ". Northeast is at top left. This figure is 
included as a separate JPEG file. 
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Fig. 3.— Selected K' (Mrk 231 and IRAS 15206+3342) and H-band (Mrk 463) data deconvolved 
with the Richardson-Lucy algorithm to show high spatial frequency structure near the nuclei. The 
indicated knots are spatially coincident with the optical knots of Surace et al. (1998), and have 
been labeled with their optical names. The one new infrared knot is labeled "irA". Northeast is at 
top left. This figure is included as a separate JPEG file. 
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optically are clearly detected in Mrk 231, Mrk 463, and IRAS 15206+3342. Figure 3 displays 
these high spatial-resolution, high surface brightness features, which have been enhanced by 
Richardson-Lucy deconvolution. These bright, compact emission regions are spatially coincident 
with the optical knots, and it seems likely that the near- infrared emission also arises from 
these knots. The many knots detected optically in IRAS 08572+3915w were not detected in 
the near-infrared, despite the near-infrared observations being deep enough to see well- defined 
tidal tails. It seems likely that this is due to a combination of heavy obscuration towards the 
infrared nucleus and the blueness of the optical knots. Even at optical wavelengths, this region's 
appearance at B and I is so different as to preclude easy registration of features. At near- infrared 
wavelengths it's appearance is that of a single source unresolved at K' ; the optical knots are too 
faint to avoid confusion with the bright nucleus. 

Perhaps most surprising is that almost no new knots have been discovered in the warm 
ULIGs. It was argued in Paper I that the optically observed star-forming knots were insufficiently 
luminous to be a significant contributor to the high bolometric luminosity in these ULIGs. 
However, one possibility is that a significantly larger population of star-forming knots exists, yet is 
embedded deeply enough so as to be undetectable in the WFPC2 data (limiting magnitude me < 
27). Assuming that most of the knots have a similar intrinsic luminosity, any knot extinguished 
by more than Ay= 4 magnitudes would have been undetectable by HST. However, the extinction 
at K'is only 1/12 that at B; knots which were extinguished just enough to be undetectable at 
B would be extinguished by only 0.3 magnitudes at K 7 . Given the K' detection limit and an 
expected K' magnitude of 20 for the brightest knots, similar knots would be seen even if they were 
obscured byAy = 25 magnitudes. This leads to the conclusion that either there is no population 
of moderately embedded star-forming knots in these galaxies or that any moderately embedded 
clusters are intrinsically less luminous than those seen optically. If any more knots exist like those 
seen optically, then they have line-of-sight extinctions greater than Ay = 25 magnitudes. The 
failure to find any large populations of spatially distributed, moderately embedded knots lends 
credence to the idea that the dust-clearing time in the vicinity of the knots is short compared to 
the current starburst lifetime. It additionally supports the finding of Paper I that the knots do 
not contribute greatly to the bolometric luminosity, since there are relatively few of them, at least 
to an optical depth of Ay= 25 magnitudes. 

An examination of the radial profiles in both bands (Figure 4) is particularly instructive. 
In nearly all cases, the ULIGs are extended on large scales (1 "~ 1 kpc) at H. The underlying 
old stellar population of the merger remnant thus seems to contribute strongly at H. The only 
notable exceptions are three systems already known to be dominated even at optical wavelengths 
by an unresolved central source: Mrk 1014, IRAS 07598+6508, and 3c273. All three are QSOs 
(Schmidt & Green 1983), and even in these cases, all three are known to possess sizable host 
galaxies (McLeod & Rieke 1994b, Boyce et al. 1996). The high luminosity of their central peaks 
is sufficient to effectively mask the low surface brightness hosts at the depth and radii illustrated 
in Figure 4. In contrast to H, at K' many of the ULIGs are dominated by a single, essentially 
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unresolved point source. This difference in compactness between H and K' is particularly evident 
in IRAS 05189-2524, IRAS 08572+3915w, IRAS 12071-0444, Pks 1345+12w, and Mrk 463e. The 
three QSOs (Mrk 1014, IRAS 07598+6508, and 3c273) show this effect to a lesser degree as they 
are already dominated by a central point source at H and K' . Of the remaining galaxies, I Zw 1, 
also a QSO (Schmidt k, Green 1983), is superimposed on an extremely large and luminous host, 
and hence has a considerable host contribution even at small radii. The Seyfert 1 galaxy Mrk 231 
is similarly known to possess a central "monster", yet it is also superimposed on an extremely 
luminous host. This leaves IRAS 15206+3342 as the only warm ULIG which is either not already 
known to be dominated by a central point-like source, or does not show increasing compactness 
at long wavelengths indicating that the majority of the bolometric luminosity is probably also 
compact. However, a single peak (the putative nucleus in Paper I) does have noticeably redder 
near-infrared colors than the other dozen or so detected knots. This may suggest that while this 
"nucleus" does not yet dominate as strongly at 2.1/xm as in the other ULIGs, it may begin to do 
so at slightly longer wavelengths. 

The physical FWHM of the PSF at K' ranges from 300 pc in Mrk 231 to nearly 1.5 kpc in 
Mrk 1014, with a median of 600 pc. In most cases the emitting region at K' (which appears to be 
spatially coincident with the optical nucleus) is at least this small. This is not surprising since 
the emission region is expected to be very small if the ULIGs harbor a dust-enshrouded QSO. 
Similarly, the ultracompact starbursts postulated by Condon et al. (1991) are also smaller than we 
can resolve in the near- infrared, and the putative nuclei identified in Paper I also have size scales 
much smaller than the near-infrared resolution limit. This result implies that the majority of the 
K'-band flux in the warm ULIGs is emitted from a region spatially distinct from the H- band 
and optical emission. This increasing compactness at long wavelengths indicates a strong infrared 
excess on small spatial scales, which we show in §3.3 is most easily explained by the presence of 
hot dust. 



3.2. Near-Infrared Luminosities 

The photometry for the total integrated galaxy light and the putative nuclei are presented in 
Table 1. The integrated galaxy magnitude represents all of the light from both the galaxy and 
the putative nuclei. It was measured using a 30" diameter aperture, corresponding to 30-100 
kpc; in all cases this was sufficiently large that the observed galaxy emission at the edge of the 
aperture was no longer visible against the background noise. The H-band luminosity of the 
putative nucleus was determined in two ways. First, detailed two-dimensional PSFs were fit to 
the galaxies; the resulting fitted PSFs were integrated in order to estimate the nuclear flux. This 
technique has several shortcomings; in particular, the fit is complicated by an ignorance of the 
underlying intrinsic galaxy profile. As an estimate, the fit was made by scaling the PSF until the 
galaxy profile reversed itself, i.e. a "hole" was dug in the nuclear regions. The other technique 
was simply to measure the total flux inside a fixed 2.5 kpc diameter aperture. Although this 



-12- 




J I I L 



12 3 4 



0.0 1.0 2.0 3.0 



i 1 1 

IR1 2071 -0444 



= I ^ I L=l L- I I I L 



i 1 1 r~- 

IR12071-0444 



12 3 4 



12 3 4 




_l I I L: 

12 3 4 




r i 

\ 

- \ 



T 1 

IR01 003-2238 



\t I I b | | | Ld 



i 1 1 r~- 

IR01003-2238~ 



12 3 



12 3 4 




~i 1 1 

IR051 89-2524. 




J I Li L 



12 3 4 



12 3 4 



i 1 1 

IR08572+391 5 




J C I L 




1 1 1 

IR08572+3915- 



_. i , tU 



0.5 1.0 1.5 2.0 



0.0 1.0 2.0 



i 1 1 r 



J I I Ld b I I I Ld 



i 1 1 r 

3C273 - 



12 3 4 



12 3 4 




1 1 1 T=l 



PKS 1 345412 



12 3 4 



J I I L: 

12 3 4 




12 3 4 



Fig. 4. — Radial profiles at H (left) and K' (right) of the warm ULIGs (eastern nuclei=filled 
circles; western nuclei=open circles) and the associated point spread functions (triangles), plotted 
as normalized surface intensity (in magnitudes) versus radius (in arcseconds). The galaxies are all 
noticeably extended at H, while they are essentially unresolved at K. Note that in terms of surface 
brightness, FWHM intensity occurs at just 0.3 mag/" below the peak. 
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Fig. 4. — (continued) 

technique ignores the underlying nuclear galaxy emission, it is fairly successful in situations where 
the central regions are dominated by a bright point source. In these cases, the surface brightness 
of the underlying galaxy is so low inside the fixed aperture compared to that of the nucleus such 
that it represents only a small fraction of the total flux. It was found that in nearly all cases, 
the two techniques agreed with each other at approximately the 10% level. Comparison of the 
photometry with that of previous studies (Sanders et al. 1988, Mcleod et al. 1994a, 1994b) shows 
that it is the same to within 15% . 

The emission at H-band is dominated by the old stellar population. The presence of double 
nuclei and well-formed tidal tails in many of these systems indicates that the merger process is not 
yet complete, and simulations indicate that they are unlikely to be more than 1-2 Gyr (Barnes 
1992, Mihos & Hernquist 1994) in age. Therefore the H-band emission is indicative of the current 
state of the stars which predated the merger itself, and the integrated H-band luminosity can be 
used to characterize the progenitors' masses by assuming a common M/Ln ratio and comparing 
them to a typical spiral galaxy. Characterization of the mass of the underlying "host" in the 
warm ULIGs is interesting for several reasons ( here we define the word "host" to mean the 
smooth component of the galaxy luminosity profile which is unlikely to be a result of clustered 
star- formation or non-stellar nuclear activity). In particular, it will allow us to determine if the 
ULIGs represent the collision of typical gas-rich spirals (Sanders et al. 1988a), or whether they are 
peculiar. For example, the ULIGs might be scaled- up versions of their less luminous counterparts, 
and simply represent the chance collision of unusually massive galaxies, in which case we would 
expect to find an unusually luminous host galaxy at H. Similarly, it may be that in order to drive 
enough material efficiently into the merger nucleus to trigger ULIG activity, the colliding galaxies 
must have very large masses. This is also important in testing any possible evolutionary scenario 
linking ULIGs and optically-selected QSOs; if QSOs are the evolved remnants of ULIGs, then 
they must have similar total old stellar masses and hence H-band luminosities. 

Table 1 presents the derived photometric quantities: the nuclear fraction of emission at H 
and K' , and the H-band luminosity. McLeod &; Rieke (1994a) estimate that for an L* galaxy 
M H =-23.9; adjusting to H =75 kms' 1 , this is M H =-24.04 . Using M v =-21.14 for an L* galaxy 
with iio=75 kms -1 and an estimate of (V-H)=2.7 for an aged stellar population with a Salpeter 
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Table 1. ULIG Total and Nuclear Photometry 



Name 


TTln, total 


TlH, nuclear 


^nuclear 
£ total 


MH.host 


m K > itota i 


, nuclear 


^nuclear 
-C-total 


I Zw 1* 


11.19 


11.69 


0.63 


-24.5 


9.94 


10.52 


0.53 


IRAS 01003-2238 


15.03 


15.80 


0.49 


-22.7 


14.30 


15.33 


0.39 


Mrk 1014 f 


12.84 


13.53 


0.53 


-25.6 


11.67 


12.28 


0.57 


IRAS 05189-2524 


11.21 


11.93 


0.52 


-24.2 


10.23 


10.44 


0.82 


IRAS 07598+6508 


13.01 








10.45 






IRAS 08572+3915W 


13.65 


15.17 


0.25 


-22.9 


12.73 


13.75 


0.39 


IRAS 12071-0444 


14.28 


15.35 


0.37 


-23.9 


13.42 


13.91 


0.64 


3c273 f 


10.88 


11.04 


0.86 


-26.1 


9.70 


9.85 


0.87 


Mrk 231 


9.91 


10.13 


0.82 


-24.4 


8.87 


8.98 


0.90 


Pks 1345+12W 


12.99 


15.36 


0.11 


-25.5 


12.62 


14.21 


0.23 


Mrk 463E 


11.19 


12.05 


0.45 


-24.8 


10.20 


10.50 


0.76 


IRAS 15206+3342 


14.14 


15.40 


0.29 


-24.1 


13.47 


14.58 


0.36 



Note. — T calibrated using photometry of Neugebauer et al. 1987. IRAS 07598+6508 of poor quality 
. Uncertainty in total magnitude is 0.05 magnitudes, except for those marked with a f, where it is 0.1 
magnitudes. Uncertainty in the nuclear luminosity is 0.1 magnitudes. 
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IMF (BC 93), we similarly derive M#=— 23.84. The range in host galaxy luminosities is ~0.4 L* 
(IRAS 01003-2238) to nearly 7 L* (3c273), with a mean of 2.4 L*. Six of the ULIGs have hosts in 
the range 1-2 L*, while three (2 QSOs and a powerful radio galaxy) are in the range 5-7 L*. This 
seems to be consistent with the idea that most of the ULIGs are the result of mergers between 
typical (i.e. L*) or perhaps slightly more luminous than normal galaxies. Previous studies have 
derived similar results: ^2.5 £g (Armus 6t si. 1990, J6ns6n ct si. 199T) sncL 2.5 ^*k' (Jgiisbii 6t 
al. 1997) for the BGS ULIGs and "Arp 220-like" objects, and ~3 L* (ranging from ~ 1-24 L*) 
for the ULIG 2 Jy sample (Murphy et al. 1996). Finally, the range and distribution in Mh for 
the "warm" ULIGs is very similar to that of the host galaxies of PG QSOs (McLeod & Rieke 
1994a, b), which typically range from 0.5-7 L^, with a mean of 1.5 L^. Additionally, McLeod & 
Rieke (1994a,b) note that the host galaxies of luminous quasars are more similar to 2 Lj^, and are 
more likely to have tidal features indicative of past galaxy-galaxy interactions. 

IRAS 08572+3915 is surprisingly underluminous («0.5 L*) at H considering that it's 
morphology is clearly that of two merging spiral galaxies (Surace et al. 1998). A possible 
explanation for this is the extremely high optical depth throughout the western galaxy component. 
Very red near-infrared colors (0.7 < H — K 1 < 1.4) within 2 kpc of the even redder nucleus 
(H-IC >2) indicate that dust obscuration may be very spatially widespread. Assuming a typical 
old stellar population color of (H— i^)=0.16, then this would imply An ~ 2-3. Dereddening this 
region by 2 magnitudes at H would increase the total galaxy brightness by 1.2 magnitudes, making 
the estimated luminosity of the underlying old stellar population much more similar to the range 
found for the other ULIGs at « 1.2 L*. We note that the derived luminosity for IRAS 08572+3915 
in K'-band is nearly 0.7 L* , further indicating that substantial extinction may obscure the old 
stars. 

This illustrates an important point: the luminosities given in Table 1 are lower limits to 
the true H-band luminosity. As illustrated by IRAS 08572+3915 above, there may still be 
considerable extinction at H due to spatially extended dust. This may affect all the ULIGs which 
have large, spatially extended regions that are redder than expected for late-type stars: notably, 
IRAS 05189-2524, IRAS 08572+3915, IRAS 12071-0444, and Mrk463. IRAS 08572+3915 is the 
most extreme example. There may also be very low surface brightness, very extended emission 
on physical scales greater than 50 kpc (e.g. in tidal tails) which we fail to detect in the near- 
infrared. This can be evaluated by comparing the spatial extent of the near-infrared data to 
previous optical data ( Sanders et al. 1988a, Surace et al. 1998). In the cases of IRAS 05189— 2524, 
IRAS 12071-0444, Mrk231, and IRAS 15206+3342 there is extended structure seen optically 
which cannot be detected in the near-infrared data. In the remainder of the galaxies, all of the 
optically detected extended structure is detected in the near-infrared, and hence this is probably 
unlikely to affect them. Finally, the techniques used to determine the nuclear emission component 
will always result in over-subtracting from the underlying galaxy. This effect is most pronounced 
when the contrast between the nuclear component and underlying galaxy core is minimal, for 
example as in IRAS 05189-2524, IRAS 08572+3915, IRAS 12071-0444, and Pks 1345+12. 
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It is interesting to note that none of these sources of error is likely to affect IRAS 01003—2238. 
This object is very blue and pointlike, does not have a high surface brightness background host, 
and has only a few knots of star formation visible optically, with no extended tidal features. The 
presence of strong Wolf-Rayet emission possibly associated with the optical star-forming knots 
would seem to preclude this being a merger so advanced as to have no visible tidal remnants. 
It is under-luminous at H (« 0.3-0.5 L*), and this, combined with the above, would seem to 
indicate that unlike the other ULIGs it is either not a merger system or it is a merger involving 
under-luminous ( and probably low-mass) galaxies. 

3.3. Optical/Near-Infrared Colors 

It is possible to examine the spectral energy distribution in the optical/near-infrared as 
characterized by the (B—I), (I—H), and (H—IC) colors by combining the new near-infrared 
data with the optical data of Paper I. This makes possible the disentanglement of the effects of 
extinction and dust emission. As has been noted previously, 2-color diagrams suffer a degeneracy 
in that many different combinations of extinction and thermal dust emission can produce similar 
results thus making it difficult to determine uniquely the underlying emission mechanism. This 
additionally complicates questions such as the age-dating of knots based on colors alone, since 
even modest amounts of reddening can result in large variations in age estimates. 

3.3.1. Modeled Colors 

The stellar colors are based on the spectral synthesis models of Bruzual & Chariot (1993; 
the updated version used is hereafter called BC95). The particular model we have chosen is an 
instantaneous starburst with a Salpeter IMF and lower and upper mass cutoffs of 0.1 and 125 
Mq. Of particular importance is that stellar (H—K') colors reach a maximum of roughly 0.16, 
which corresponds to an extremely old age (> 10 9 yrs) for a stellar ensemble. This is very similar 
to the observed colors of typical galaxies (Aaronson 1977). 

The QSO colors are based on a model of a typical QSO optical/near- infrared SED. The 
continuum emission is modeled with a two-part broken power law derived from Neugebauer et al. 
(1987) for the Bright Quasar Survey (BQS): 

F A oc A~( 2+a ) where a = —1.4; A > 1/xm 

a = —0.2; A < 1/xm 

Broad line emission based on the line widths, equivalent widths, and line strengths determined 
by Baldwin (1975), Davidson & Netzer (1979), Wills et al. (1985) and Osterbrock (1989) was 
also modeled. The derived QSO colors as a function of redshift are given in Table 2. The colors 



(3) 
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derived from this model are in good agreement with the observational results of the Elvis (1994: 
UVSX) sample, whose average colors are (fl-J)=l.l±0.3, (I- fl)=1.3±0.3, and (H-K' )=0.78±0.2 
at a median redshift of 0.15. The presence of strong line emission is the cause of most of the 
dependency of the color on redshift. In particular, Ha, H/3, and Mgll are shifted through the B 
and I filters in the redshift range < z < 0.16. The considerable scatter in QSO colors is likely to 
be intrinsic. Neugebauer (1987) noted a considerable scatter in the continuum power-law indices; 
these alone can easily account for the observed color variances. 

The effects of line-of-sight extinction, i.e. a simple foreground dust screen, were derived 
from Rieke & Lebofsky (1985). The reddening law was linearly interpolated to the rest frame 
wavelength of the galaxy in order to account for redshift. 

The spectrum of hot dust with a characteristic temperature, T, and a A" 2 emissivity law was 
modeled according to the prescription given by Aaronson (1977): 

F x oc \~ 2 B X oc \- 7 (E^r - l)- 1 (4) 

Modeled dust temperatures were in the range 600-1000 K. The effects of cooler thermal dust 
emission at long wavelengths on color are degenerate since the dust emission is so red that it 
only increases the K'-band flux, and hence the (H— IC) color. Only at temperatures above 600 
K can dust emission contribute appreciably at wavelengths shorter than 2 fim. 1000 K is near 
the realistic ceiling for dust temperatures; at temperatures above ~1500 K, the dust evaporates 
(Spitzer 1978). In any case, such a hot dust model does not fit the data well as even 1000 K is 
clearly hotter than needed. 

Free- free emission was modeled according to (Spitzer 1978) 

F x oc A~VaW (5) 
assuming a 20,000 K electron temperature. The A -2 term dominates this expression; the slowly 



Table 2. Derived QSO 
optical/near-infrared Colors 



z 


(B-I) 


(I-H) 


(H-K') 


0.00 


0.54 


1.75 


0.78 


0.05 


0.61 


1.68 


0.78 


0.10 


0.75 


1.50 


0.78 


0.15 


0.93 


1.20 


0.78 
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varying Gaunt factor is regarded as constant. Similarly, at longer wavelengths the exponential 
term is near unity, and only the short wavelength colors are affected appreciably by it. However, 
these same short wavelength colors will be shown to be poor at discriminating the underlying 
emission mechanism. 

Finally, the case of mixed stars and dust was considered. The spatial extent of the knots 
is such that they must be very large associations of many stars. It is both possible and likely 
that there are local, unresolved regions of high extinction inside the knots. An example of this 
might be a small region of star-formation containing several embedded protostars. The result of 
the presence of such unresolved regions inside the knots is to increase the thermal near-infrared 
component of the SED relative to that of a naked starburst, in a manner similar to that of the hot 
dust modeled above. Note that this is different from the effects of uniform line-of-sight extinction 
as modeled above. 

We have modeled the effect of mixed stars and dust by assuming an underlying power- law 
distribution for the unreddened luminosity La v obscured by a given extinction Ay'. 

L Av = L A=0 (Ay + l) a (6) 

This is an "onion-skin" model; it peels away each optical depth and computes the total luminosity 
in each filter from all emission at that optical depth. The colors were derived by computing a 
numerical model of the observed luminosity at each optical depth, and then summing over all 
extinctions from Ay=0 to Ay=A mSuX . The resulting colors are displayed in Figure 5. The a=0 case 
corresponds to dust mixed uniformly with the emitting source, and is similar to that of Aaronson 
(1977) and Carico et al. (1990). The {B—T) colors rapidly converge on a fixed value determined 
by a, since the effective optical depth rises very rapidly at short wavelengths and hence stars at 
Ay >10 magnitudes cannot contribute more than a tiny fraction to the observed short wavelength 
emission. At Ay > 50 magnitudes this becomes true even at K' , and thus each model converges 
at high optical depth on a color dictated by a. This behaviour can be most easily described by 
the statement that the colors are fixed by the total stellar luminosity at optical depths less than 
one for any given wavelength, i.e. only the relatively unobscured stars in the outer regions of 
the dust/star mixture contribute noticeably to what is observed. This is similar to the result of 
Witt et al. (1992) that any extended stellar distribution mixed with a scattering and absorbing 
interstellar medium reaches a maximum reddening value at sufficiently high optical depths. While 
models with large amounts of stars at high optical depth (e.g. a=2) have fairly large {H—K 1 ) 
colors, they are accompanied by a similar increase in (B—T) and (I—H) which is not entirely 
dissimilar to the effects of reddening. In the reddening-orthogonal color basis discussed below, 
nearly all four models lie in essentially the same color location, and hence only the a=0 case will 
be illustrated. We have assumed that the luminosity source has the same colors as an intermediate 
age starburst (B—I, I—H, H—K'= 0.62, 0.97, 0.2), and that the dust obeys the same reddening 
law discussed above. The former assumption is almost certainly unrealistic, since it ignores the 
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thermal dust emission associated with protostars found inside star-forming regions. However, this 
will be discussed using an additional observational model presented in §3.4. 

The treatment described above ignores the effects of scattering, which may be quite important 
for an extended distribution of stars and dust. Witt et al. (1992) present models for various 
spherical distributions of dust and stars, with a full treatment of scattering. Of particular interest 
here are those models for "dusty" galaxies in which dust and stars are uniformly mixed. This 
corresponds to the model generated above, and is the most likely to represent the situation within 
the starburst knots (other models which represent dusty starburst cores or obscured AGN are 
more similar to the galaxies taken as a whole, rather than as small-scale individual knots). The 
addition of scattering leads to an increase in the fraction of emergent light, particularly at short 
wavelengths. The evolution of optical/near-infrared colors for both models is similar under the 
assumption of an intermediate age burst as described above. With a total extinction of Ay = 16.2 
magnitudes, the non-scattering model predicts (B—I, I—H, H—K f ) = (1.63, 1.96, 0.54), while the 
scattering model predicts (B-I, I-H, H-K f ) = (1.46, 1.60, 0.51). 

A more important change, however, is that scattering reduces the effective optical depth by 
increasing the fraction of light originally emitted by the stars which actually emerges from the 
model. This is particularly valuable because the mixed stars and dust model without scattering 
developed above cannot easily constrain the total bolometric luminosity of the stellar ensemble; it 
is relatively easy to hide luminosity at high optical depth where it cannot be seen in comparison to 
the stars at lower optical depth. One can therefore easily increase the "geometric" correction (the 
difference in magnitude between the true luminosity and that observed) at any given wavelength 
by simply specifying a sufficiently large total optical depth for the star /dust mixture. In the 
uniformly mixed stars and dust model without scattering and with a total optical depth of Ay 
= 50 mag, for example, the geometric correction to convert between the actual total stellar 
luminosity and that observed is nearly 4.8 magnitudes at B and 3.5 at K. For an ensemble with Ay 
= 16 magnitudes these numbers drop to 3.1 and 0.8, respectively. For the "dusty" galaxy model 
with scattering, however, even if the total optical depth at B is as high as r=20 (corresponding 
to the Ay =16 magnitudes case described above), the geometric correction is only 2.4, and at 
K is 0.7. This effect is even more pronounced for other dust geometries; at Ay =16 magnitudes 
the geometric correction at B is 0.4 for stars with dust lying interior to the stars and 1.8 for a 
"starburst" system with a centrally concentrated stellar distribution and dust mixed with the 
stars. Witt & Gordon (1996) showed that this effect becomes even more pronounced when the 
scattering medium is assumed to be clumpy instead of homogeneous, in which case the effective 
optical depth may be decreased by additional factors of 2-4. The addition of the scattered light 
to the total observed light therefore makes the observed luminosity more indicative of the actual 
total luminosity. 

Figures 6 and 7 show the three-dimensional space occupied by the three colors (B— I), 
(I—H), and (H—IC ). The resulting cube has been rotated by roughly 90° in order to illustrate the 
orthogonal nature of the extinction and dust emission curves. Note that stellar and QSO nuclear 
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Fig. 5. — {I—H) and (H—K* ) colors for models of mixed stars and local dust extinction, with the 
luminosity at each optical depth characterized by a power-law of index a. The marked intervals 
are the maximum optical depth in units of Ay = 10 magnitudes. The (B—I) colors for a given 
model become fixed rapidly since only stars at low optical depth contribute significantly to the 
short wavelength emission. This is true even at K'-band at high optical depth, hence the colors 
converge on a model-dependent value. The open circle indicates colors for mixed star and dust 
models including scattering with Ay = 16 magnitudes, although (B— I) is lower (1.46 for evenly 
mixed stars and dust, and 0.77 for centrally concentratated stars and dust). The same point in the 
non-scattering models is marked with a small filled circle. 
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emission are well separated in the three-color representation. The "dusty" galaxy scattering 
model occupies nearly the same color space as the free-free emission model. The (B—T) colors 
of QSOs are similar to the expected colors of moderate age (10-100 Myr) starburst clusters; the 
largest distinguishing characteristic between the two is that QSO nuclei are more luminous at long 
wavelengths. As a result, the primary advantage of the short wavelength {B—T) data is to allow 
the disentanglement of extinction effects. 

The following discussion will concentrate primarily on Figure 7, which has been rotated such 
that the color basis projected onto the page is orthogonal to extinction (i.e. line-of- sight extinction 
does not change the apparent position of the points). Since extinction is almost certainly present 
to a varying degree throughout these galaxies, it is important to eliminate as much as possible 
extinction effects. To distinguish between stellar emission and QSO activity, any processes capable 
of producing similar optical/near-infrared colors must be considered. In the following discussion, 
"reddening" will be used to mean processes which move points upward in the projected basis of 
Figure 7, primarily giving them larger values of (H—IC), as opposed to the reddening caused by 
extinction. There are few processes which can make the QSO emission bluer. Large amounts of 
free- free emission (open circles, the last of which is pure free- free emission) can decrease (H—K 1 ) 
by only as much as 0.2. Likewise, free-free emission can redden the stellar colors only slightly. 
The most important effects are clearly that of hot dust and dust mixed with stars. The overall 
effect of mixed dust and stars is to redden the basic stellar colors. As discussed above, however, 
this reddening saturates when sufficiently high optical depths are reached. Furthermore, models 
of distributed stars and dust with both scattering and absorption indicate that it is difficult to 
redden starlight even with very high total optical depths (Witt et al. 1992,1996), and hence it is 
difficult to mimic the very red near-infrared QSO colors. 

Since the QSO and stellar colors differ primarily at long wavelengths, the most serious 
problem occurs with reddening processes that affect only K' and longer wavelength bands, such 
as thermal dust emission. If 800 K dust contributes as much as 55-60% of the emission at K' , 
then the resulting stellar /hot dust ensemble will have colors nearly identical to that of optically 
selected QSOs. Thankfully, some reprieve is noted in that the extremely red colors of hot dust 
can also redden the QSO nuclei, and hence any hot dust component around a QSO moves it even 
further away from the basic starburst colors. Unfortunately, the very red dust colors also imply 
that the stellar-hot dust emission will converge on the same colors as QSO-hot dust emission as 
the percentage of hot dust increases. For stellar colors, this effect is noticeable when hot dust 
accounts for > 70% of the emission at K' . Carico et al. (1990) found that for far-infrared bright 
and luminous (/6o > 5 Jy; > 10 11 Lq) galaxies which do not show signs of Seyfert activity, 
typically {H-K 1 )= 0.40 (ranging from 0.16 to 0.57). Since these galaxies are believed to be powered 
predominantly by star formation, this represents the typical near-infrared colors observed for the 
most powerful of clearly identifiable starbursts. 

Because the detailed SEDs of the high spatial frequency structure in the warm ULIGs are 
unknown, particularly at short rest wavelengths, it is not possible to apply K- corrections to 
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Fig. 6. — (B—T), (I—H), (H—K f ) color cube illustrating the colors of the following emission sources: 
(1) optical QSO nuclear colors indicated by the large, dotted circle and (2) the evolutionary track 
of an instantaneous starburst aging from to 15 Gyrs, indicated by the open circles. It also shows 
the effects of the following physical processes on the colors of stars and QSOs: (1) line-of-sight dust 
extinction in units of Ay = 1 magnitude (open squares), (2) free- free emission with an electron 
temperature of 20,000 K in increments of 20% of the total flux at K' (open joined circles), (3) 
emission from 800 K dust in increments of 10% of the total flux at K' (two sets of filled, joined 
circles representing hot dust combined with either a QSO or young stellar SED), and (4) emission 
from uniformly mixed stars and dust, in units of Ay =10, 30, and 50 magnitudes (x). Note that 
the line-of-sight dust extinction and thermal dust reddening curves are nearly orthogonal. 
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Fig. 7. — Same as Figure 6, but rotated such that the projected color basis is orthogonal to the 
line-of-sight dust extinction, i.e. line-of-sight extinction does not affect the position of points in this 
graph. Note the separation between QSO and stellar colors. Also, nearly all the warm ULIG nuclear 
colors are very similar to that of the infrared-loud QSOs, with the exception of IRAS 01003— 2238. 



-24- 



them. However, since the models of the emission processes, starbursts, and QSOs do have detailed 
SEDs, inverse K-corrections can be made to their rest-frame colors at the redshift of our targets. 
This is done by convolving the synthetic spectra with the known detector and filter bandpasses. 
The magnitude zeropoint calibration of the filters is derived using the Kurucz model spectrum of 
Vega (BC95). For brevity, Figures 6 and 7 are calibrated to the median ULIG redshift, z=0.12. 
The K-corrections affect the modeled stellar colors by 6{B- 1,1- H,H- IC )=(0.05,0.07,0.18) for a 
young (10 Myr) starburst, and by (0.40,0.11,-0.08) for an old (2 Gyr) population. The effects of 
K-corrections can therefore be quite large depending on the modeled population, and hence the 
representation of stellar colors can be considered to have a sizable uncertainty attached. Larger 
redshift values also shift the dust emission curves to a more vertical orientation, as the rest frame 
filter bandpasses become bluer. 



3.3.2. Nuclear Colors 

Also shown in Figures 6 and 7 are the colors of the ULIG nuclear regions (open triangles). 
IRAS 01003—2238 seems to be unique in that it is most consistent with a young (10 7 yrs) stellar 
population mixed with local dust extinction with an optical depth of Ay magnitudes. It could also 
be a combination of hot young stars with 20% hot dust at K' and Ay = 1 magnitude of extinction. 
IRAS 15206+3342 has an SED similar to that of a QSO reddened by Ay = 1 magnitude, but it 
is also similar to a young stellar population with a very large thermal dust component and Ay = 
1-2 magnitudes of line-of-sight extinction. Given the large range of QSO colors it is not possible 
to differentiate between these two possibilities based on the optical/near- infrared colors alone. 

The remaining 9 galaxies all have SEDs characterized by large K' excesses. The three 
optically identified QSOs (IZwl, Mrkl014, and 3c273) may be taken as representative of a 
subclass of "infrared loud" QSOs. Mrk231, the most QSO-like of the ULIGs, has an SED nearly 
identical to that of the infrared-loud QSOs. Similarly, IRAS 05189-2524, IRAS 08572+3915w, 
and IRAS 12071—0444 also have colors essentially identical to that of the infrared-loud QSOs seen 
through Ay ~ 2 magnitudes of extinction. PKS 1345+12w has a smaller near- infrared excess 
than the infrared-loud QSOs, but has the same colors as a typical QSO (characterized by our 
synthetic colors) seen through Ay = 2 magnitudes of extinction. Mrk 463e has a larger excess at 
both H & K'than the others. This can be explained by emission from 800-900 K thermal dust 
which contributes 20-30% more of the K- band flux than exists in the infrared-loud QSOs. 

Compared to the synthetic optically-selected QSO colors, both the ULIGs and the subsample 
of three infrared-loud QSOs have a noticeable K'-band excess. This excess is 2-3o" greater 
than that expected considering the normal range of QSO colors (<7fl_x'=0.2; Elvis et al. 1994). 
However, it is likely that previous near- infrared and optical surveys have underestimated the 
(H—K t ) color of the QSO nuclei. Given that normal galaxy colors (H—K 1 <0.2; Aaronson 1977) 
are almost always bluer than QSO nuclei (H—IC>0.8), studies that used large apertures to 
measure the QSO nuclear contribution are contaminated by underlying host galaxy light and 
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therefore derive colors bluer than the actual nuclei. Also, several authors (Neugebauer et al. 1987, 
Sanders et al. 1989) have suggested that a hot dust component already contributes significantly to 
the near-infrared emission in optically selected QSOs. Given the presence of a strong far-infrared 
dust emission component and significant molecular gas content (see Sanders 1991), it would not 
be surprising if a larger reservoir of dust has also given rise to enhanced hot dust emission. 

3.3.3. Knot Colors 

In Table 3, we present photometry for the additional compact sources observed in the 
warm ULIGs. This includes both the galaxy "nuclei" which do not appear to have AGN- like 
characteristics (e.g. IRAS 08572+3915e) and are likely to be the true stellar cores of galaxies, 
and the clumps of near-infrared emission which we associate with the compact star-forming 
"knots" observed in Paper I. The galaxy nuclei were measured using aperture photometry on 
the undeconvolved images, and have an uncertainty of 0.1 magnitudes. The sky background 
was determined by using annular regions outside the aperture which appeared free of high 
frequency structure. The knot magnitudes were measured from the R-L deconvolved images 
in order to minimize confusion. Uncertainties in the photometry are introduced by the non- 
linear, non-invertible nature of the deconvolution algorithm; these were checked by comparison 
with aperture photometry of the raw images. The knot magnitudes are uncertain by 0.15- 0.2 
magnitudes. 

The colors of the starburst "knots" found in Paper I were examined and initially revealed 
anomalous colors: while the knots had (B—I) and {H—K 1 ) colors that were consistent with stars, 
the (I—H) colors were far larger (typically (I— H)> 2.5) than could be explained by any plausible 
combination of stellar colors and reddening effects. A careful examination of the data revealed 
that this was most likely due to aperture effects. The resolution of the optical data is typically 
5-10 times higher than that of the infrared data. As a result, the apertures used to measure the 
optical data are much smaller. It is likely that the larger apertures used for the near-infrared data 
also include some additional underlying emission that may not be related to the optical knots, 
hence the apparent increase in luminosity between the optical and near-infrared data. However, as 
noted previously, the largest differences in emission and obscuration mechanisms lie in the optical 
(B—T) and near-infrared (H—K' ) colors. The purely optical and purely infrared knot colors should 
be fairly certain since data of different resolutions are not mixed. Therefore, the knot colors in the 
2-color (B— I) vs. (H—K) plane are presented in Figure 8. 

The colors of the galaxy nuclei appear very similar to that of stars with mild reddening. 
Pks 1345+12e has colors similar to that of an old stellar population with Ay=l-2 magnitudes. 
IRAS 08572+3915e and Mrk 463w appear similar to a somewhat younger stellar population (1-2 
Gyrs), but have a slight near-infrared excess. In the case of Mrk 463w this may not be surprising 
as this galaxy core shows both Seyfert 2 activity and numerous knots of star formation near it, 
indicating that there may be considerable new star formation in it, resulting in a thermal dust 
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Table 3. Additonal Structure 



Name ra H m K ' 



IRAS 08572+3915 
^(east) 1 " 16.2 16.0 
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Pks 1345+12 


6(east) t 


15.8 


15.5 


Mrk 463 


19 


17.3 


16.8 


20(west) t 


14.3 


13.8 


21 


16.9 


16.7 


28/29/30 


18.0 


17.6 


IRAS 15206+3342 


1 


17.6 


17.4 


2 


18.0 


17.0 


3/4/5 


17.6 


17.1 


7/8 


18.1 


17.4 


gt 


16.9 


16.2 


10 


18.8 


18.1 


n 


17.9 


17.3 


12/13 


17.5 


17.2 


14 


19.1 


18.7 


15/16 


<20.6 


<19.8 


17 


18.1 


17.5 


18 


<20.4 


<19.2 


19 


18.2 


18.0 



Note. - ' denotes probable second (non-active) 
nucleus in double-nuclei systems. Names taken from 
Suracc et al. (1998) and should be referred to in the 
form (galaxy name):SSVVM(number). Names with 
multiple numbers indicate the given magnitude is an 
integrated value for the named knots. 
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Fig. 8. — The (B—T),(H—K f ) color plane, with the same reddening mechanisms, stellar and QSO 
colors as in Figure 6. The star-forming knots are marked with stars and the stellar galaxy nuclei 
are labeled. 
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contribution. There is a surprisingly wide range in (H—K 1 ) colors for the putative starburst knots 
themselves. For comparison, we consider larger samples of luminous infrared galaxies (LIGs; 
10 11 < L 1T < 10 12 Lq; Carico 1988, 1990). These are the most far-infrared luminous galaxies 
in which it is fairly certain that star formation is the primary contributor to the far- infrared 
luminosity, and as such represent the most extreme known starburst galaxies. The distribution in 
knot colors is very similar to that of the LIG nuclear regions, (H—K)= 0.25-0.6, not counting 
the few LIGs known to be Seyferts and which are therefore likely to have a strong non-stellar 
component. 

Carico et al. (1990) found that the "disk" (H—K > ) colors of luminous IRAS galaxies were 
typically (H—K )=0.25 (ranging from 0.08 to 0.37). It is unclear what stellar components were 
integrated into this "disk" light — clearly there must be some contamination by star-forming 
regions, which probably accounts for the difference between these colors and the value of 
H—K =0.16 found by Aaronson (1977). The colors of the underlying host galaxies of the warm 
ULIGs were examined by subtracting the contribution of the putative nucleus as well as any 
additional contribution from identified point-like sources from the integrated galaxy flux. In 
this case the contribution from the putative nucleus was taken as the flux inside the central 2.5 
kpc. This should ensure that all the nuclear contamination of the underlying galaxy colors is 
removed and, although this process also eliminates some of the host galaxy's flux, should leave 
the colors unchanged. Most of the colors so derived for the warm ULIGs appear similar to 
those found by Carico et al. (1990), ranging from H-K =0.09 in IRAS 05189-2524 to nearly 
H-K =0.7 in IRAS 08572+3915 (the latter probably being strongly affected by dust) with a 
median of H—K=0A. However, a few of the systems seem to show peculiarly red colors, namely 
the optically selected QSOs. This is probably the result of residual nuclear contamination caused 
by the extended wings of the bright quasar nuclei. 

3.4. Contribution of putative "nuclei" to the bolometric luminosity 

ULIGs are defined by their extremely high infrared luminosities (L\ x >10 12 Xq, equivalent 
to the bolometric luminosities of optically selected QSOs). If the source of the majority of this 
luminosity were an AGN, then it is likely that after expulsion of the enshrouding dust the AGN 
would appear similar to an optical QSO. Although the broad-band optical/near- infrared colors of 
the warm ULIG nuclei (as defined in §3.1, the dominant compact emission sources at K' spatially 
coincident with previously identified optical nuclei) appear very similar to that of QSOs, many 
of the nuclei are less luminous at these wavelengths than QSOs. Three of our 12 ULIGs (the 
optically identified QSOs IZwl, Mrkl014, and 3c273) already have a substantial fraction of their 
Lbol hi a QSO-like SED at short wavelengths (i.e. "the blue bump"). The nucleus of the broad- line 
object IRAS 07598+6508 also meets the minimum 5-band luminosity criterion for QSOs (Surace, 
unpublished). Dereddening of the observed optical emission by the optically estimated extinction 
(typically Ay = 2 magnitudes, Surace et al. 1998) is insufficient to raise the estimated optical 
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luminosity of the remaining eight AGN to the canonical Mb = —22.1 luminosity criterion for 
optical QSOs except in the case of the nucleus of Mrk 231 (dereddened Mb = —23 ). Dereddened 
B-band luminosities of the remaining seven range from Mb=— 15.5 for IRAS 08572+3915 to Mb 
= -22 for IRAS 15206+3342. The median value is M B = -20.2, which is typical of Seyferts 
(Osterbrock 1989). 

Due to the much lower extinction in the near-infrared, it is possible that the nuclei in many 
cases are seen directly in the near-infrared (which is supported by their essentially unresolved 
near-infrared morphology) and hence the observed K'-band luminosity may be more indicative 
of their true luminosity, although it is also sensitive to hot dust emission, as in §3.2. Using 
the observed colors for IZwl, Mrk 1014, and 3c273 (the infrared-loud optically-selected QSOs), 
Mb=— 22.1 is equivalent to a QSO K 7 magnitude of Mk'=— 25.8; using the synthetic optical 
QSO colors it is M K '=-25.2 at the median redshift z = 0.12. IRAS 05189-2524 and Mrk 463e 
(which are also broad-line objects) now qualify as QSOs based on their K' luminosity. The 
remaining five of the warm ULIGs still fall short of QSO near-infrared luminosities; their infrared 
luminosities range from M K '=-23 (IRAS 01003-2238) to M K '=-24.7 (IRAS 12071-0444), with 
a mean of Mk'=— 23.8. Therefore, while in 58% (7/12) of the warm ULIGs the dereddened 
near-infrared nuclear luminosity is similar to an optically selected QSO and is probably the source 
of the majority of the bolometric luminosity, if we assume only a simple foreground reddening 
screen then the remaining five warm ULIG nuclei appear similar to lower luminosity Seyferts 
and fall an order of magnitude short of QSO luminosities. It is notable, however, that in these 
underluminous systems both the nuclei and the star-forming knots are insufficiently luminous at 
optical/near- infrared wavelengths (even after dereddening) to account for the high bolometric 
luminosity. Either a large fraction of the total luminosity must arise in a source or sources so 
heavily extinguished as to escape detection at both optical and near-infrared wavelengths, or 
alternatively the "knots" /nuclei we have already detected must have intrinsic SEDs with an 
unusually strong far-infrared component in comparison to the putative nuclei in the seven ULIGs 
which have dereddened optical/near- infrared luminosities consistent with QSOs. 

Until now we have assumed that both the optical and near-infrared emission reached the 
observer through a common path. The apparent optical and near-infrared underluminosity in five 
of the warm ULIGs relative to QSOs may be understood as a putative AGN being seen through a 
combination of patchy extinction and scattering. The nucleus might in some cases be seen directly 
in the near-infrared due to the lower near- infrared extinction, but primarily through scattered 
light in the optical (which may then be additionally reddened by extinction along the scattered 
path). In more extreme cases the line-of-sight column density may be so high that the AGN is 
seen only in scattered light even at near-infrared wavelengths. Since the extinction estimates are 
most heavily influenced by the short wavelength colors (where extinction effects are greatest) the 
extinction estimates would be dominated by that along the low-extinction scattered paths, leading 
to estimated extinctions that are much too low. Attempts to derive true luminosities based on the 
optical light would then be futile if they are dominated by scattered light which is dependent on 
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the scattering geometry. 

Polarization studies (Young et al. 1996) which included 6 of our warm ULIGs have actually 
found evidence for patchy obscuration in the vicinity of the scattering region. In particular, 
their conclusions for IRAS 05189— 2524 agree with the conclusion presented in §3.1 based on 
morphology that the emitting region is seen directly in the near-infrared, but with an optical 
scattered component with mild optical extinction and a sizable contribution from starlight, which 
is probably a result of the circumnuclear starburst surrounding the AGN (Paper I) . A similar result 
is found for Mrk 463e, whose optical/near-infrared morphology shows that the optical emission is 
dominated by a spatially distinct region to the north of the near-infrared nucleus. Recent HST 
results indicate that this emission is that of a highly polarized Seyfert 1 "hidden" nucleus (Tremonti 
et al. 1997). This agrees with the polarized SED found by Young et al. (1996) who concluded 
that the line-of-sight extinction to this nucleus was Ay=17 magnitudes; this high extinction 
agrees with our observed morphology. The steeply rising polarization in the near-infrared of 
IRAS 08572+3915w (Young et al. 1996) may indicate that it is so deeply buried that almost 
none of the optical emission arises from the AGN, and that even at near-infrared wavelengths 
the emission is predominantly scattered and not line-of-sight. The case of IRAS 08572+3915w 
reinforces the point that morphologically the galaxies with the least luminous nuclei also tend to 
have more complex optical morphology, suggestive of patchy scattering. Additionally, recent high 
spatial resolution observations of nearby Seyfert galaxies have indicated the importance of patchy 
extinction in differentiating Seyfert classes (Malkan et al. 1998). Given that the far- infrared 
emission in the warm ULIGs is 50-90% of the bolometric luminosity, then this implies that the 
covering factor of a putative QSO must be of the same order. Patchy extinction must therefore 
be present; only in those cases where the nuclei are seen along a direct line of sight would the 
unreddened luminosities be consistent with QSOs. 

In terms of the previous color analysis which assumed a simple single-path model for the 
optical/near- infrared emission, correcting for the effects of patchy extinction and scattering would 
be to drive the nuclear colors even further from those of stars, thus further eliminating some of the 
degeneracy in the color diagram. The apparent enhancement at shorter wavelengths due to the 
addition of scattered light would result in the derivation of colors that are too blue — correction 
for this further reddens the colors. While most polarization studies detect only small fractional 
polarizations (< 5%), they also use very large apertures (5-8"), which encompass large fractions of 
more distant galaxies like the warm ULIGs studied here (Young et al. 1996). The beam dilution 
can be estimated by comparing the B-band emission of the putative nuclei in the HST/WFPC 
data to the total flux encompassed in a 5" aperture. If the polarized light is actually associated 
with the high surface brightness features observed with HST then the beam dilution correction 
may be anywhere from 5x in Mrk 463e to as high as 50 x in IRAS 08572+3915, with 10-15 x 
being more typical. In this case, a significant fraction of the observed optical nuclear emission may 
be scattered, affecting the analysis as mentioned previously. 
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3.5. Contribution of putative star-forming "knots" to the bolometric luminosity 

In some of the objects in our sample of warm ULIGs the starburst knots (compact emission 
regions not identified as "nuclei") may contribute as much as 25% of the bolometric luminosity. 
Since the ULIG knots appear to be powerful bursts of star-formation similar to the circumnuclear 
starbursts found in the LIGs, they may have similar SEDs. The relationship between K 7 and 
far-infrared luminosity in the LIGs can then be used to estimate a far-infrared luminosity of the 
ULIG knots based on their observed near- infrared luminosity. From Carico et al. (1990) the 
empirical relation 

M W - 6.45 

logL; r = 7 

is derived for the subset of LIGs without any evidence for a strong AGN component; this relates 
the "nuclear starburst" near- infrared luminosity to L ir . For those warm ULIGs where it is possible 
to observe bright near- infrared knots, the extrapolated total far-infrared knot luminosity ranges 
from L h = 10 10 - 27 L Q in Mrk 231 to as high as log L lT = 10 116 L Q in IRAS 15206+3342. Previous 
estimates of the knot bolometric luminosities (e.g. Surace et al. 1998) were appreciably lower 
(0.1-2.5 x 10 10 Lq) since they were based on their B magnitudes and a simple foreground 
reddening model applied to a purely stellar (no dust) starburst. The new near-infrared observations 
reveal a near- infrared excess in the knots similar to other infrared-luminous starbursts which 
results in higher estimated bolometric luminosities. Nevertheless, the net result is that the knot 
contribution to Lboi is much less than that of the putative nuclei, which are thus likely to be the 
dominant energy source in the warm ULIGs. 



4. Conclusions 

We have presented high spatial resolution images of a complete sample of "warm" 
ultraluminous galaxies. From these images we conclude the following: 

1) The warm ULIGs have a near-infrared structure very similar to that seen optically, with 
tidal features and some clumpy emission coincident with the star-forming knots seen optically. 
Most of these knots have low line-of sight extinctions. We have failed to detect significant numbers 
of new near-infrared knots indicating that any additional knots are either extinguished by greater 
than ^4y=25 magnitudes, are intrinsically less luminous, or are distributed differently from the 
optical knots. 

2) The starburst knots have a wide range of {H—K 1 ) infrared excesses. These excesses span 
the range seen in the less luminous LIGs. Adopting the SED of the LIGs as representative of the 
most luminous star-forming regions, then the knots seen in the warm ULIGs have total, estimated 
far-infrared luminosities in the range of 10 9 - 10 11 ' 5 Lq. This is similar to the range found in 
non-active infrared galaxies. 
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3) A single "knot" of emission, spatially coincident with an optically identified "nucleus", 
increasingly dominates the emission at long wavelengths. This putative nucleus is generally 
unresolved at K 7 . 

4) The colors of these putative nuclei are nearly identical to that of optically selected QSOs 
with an added hot (~800 K) dust component contributing 10-30% of the K' emission seen through 
Ay=l-2 magnitudes of line-of-sight extinction. These same colors are identical to "infrared loud" 
quasars. They could, however, also be produced by a powerful starburst 10-100 Myrs in age with 
70% or more of the flux at K' originating in hot dust. This would be a greater percentage than 
that observed in other, less luminous starburst galaxies. 

5) Seven of the sample of 12 warm ULIGs have dereddened "nuclear" optical and near- 
infrared luminosities consistent with that of optically-selected QSOs. In the remaining five cases 
the dereddened luminosities of the putative nuclei are similar to that of Seyferts. If the underlying 
emission source is a QSO, then this could be explained by the optical emission being seen primarily 
through scattered light along less extinguished lines of sight. 

6) The powerful Wolf-Rayet galaxy IRAS 01003—2238 appears unique among the warm 
ULIGs: it is pointlike, has no obvious tidal features or extended host, and is under- luminous at 
both optical/near-infrared wavelengths. 

7) The dereddened "nuclear" luminosities of the sample of warm ULIGs are sufficiently high 
that in 58% (7/12) of the cases they could provide most of the bolometric luminosity. In all cases 
except IRAS 15206+3342 their contribution to the bolometric luminosity is likely to greatly exceed 
that of the more widespread star- formation detected so far. 
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systems. We also thank Lee Armus and B.T. Soifer for helping with the early groundwork for this 
project. Finally, we thank Bill Vacca, Alan Stockton, Steve Willner, and an anonymous referee 
whose useful comments helped strengthen both the material and presentation of this paper. J.A.S. 
and D.B.S. were supported in part by NASA grant NAG5-3370. 
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A. Notes on Individual Objects 

I Zw 1 — a Palomar-Green (PG) QSO. Both spiral arms are clearly detected in the 
near-infrared. The QSO nucleus is situated on an extremely bright elongated stellar galaxy core. 

IRAS 01003-2238 — A powerful Wolf-Rayet galaxy (Armus et al. 1988) with several small 
knots of star-formation seen optically. The near-infrared data reveals no features at all. 

Mrkl014 — a PG QSO observed optically to have several knots and linear features very 
close to the nucleus (Paper I), and a very large tidal arm (MacKenty & Stockton 1984). The 
single-armed tidal feature to the east of the optical nucleus is clearly detected in the near-infrared. 

IRAS 05189— 2524 — the extended structure of the host galaxy is clearly detected at H. 
Under careful scrutiny, the eastern and southern loops can also be detected. Young et al. (1996) 
report broad polarized Ha and polarization evidence that the nucleus is seen directly in the 
near-infrared, confirmed by our detection of an apparently unresolved nucleus. This contrasts with 
the optical morphology presented in Paper I which showed a double peaked feature that appeared 
most consistent with a single object bisected by a dust lane. 

IRAS 07598+6508 — our data for this object is poor. However, we do not readily detect the 
tidal structure and star-forming knots seenoptically to the south and east of the nucleus (Boyce et 
al. 1996), probably indicating very blue colors. 

IRAS 08572+3915 — in the near-infrared, both nuclei can be clearly seen, as can the two 
tidal tails extending to the north and east. The western nucleus appears similar to a point 
superimposed on an extremely bright, extended background. The eastern nucleus is oval in shape 
and is very similar both in infrared and optical appearance. 

Mrk 231 — this system is dominated in the near-infrared by the active nucleus seen optically. 
Additionally, the optical knots in the southern "horseshoe" are detected, and are even more distinct 
in the deconvolved images. Additionally, we seem to have detected at least one near-infrared knot 
not seen optically "irA". 

Pks 1345+12 — the extended galaxy is clearly seen and is similar in appearance to previous 
optical images, and has no fine structure. The underlying host galaxy is very luminous and thus 
the nuclei contribute less than half to the total near-infrared luminosity. 

Mrk 463 — this system is dominated at K'by the eastern nucleus. The peak of the emission is 
located at knot 15 of Surace et al. (1998). This seems to imply that this knot is the true nucleus, 
and that the optically bright spot to the north is an excitation feature. This is probably the source 
of the polarized broad lines reported by Tremonti et al. (1997). The extended tidal remnants are 
detected, as are the star forming knots in these tails (primarily to the north and east). 

IRAS 15206+3342 — as noted in Surace et al. (1998), this system possesses 4 extremely 
bright optical knots, as well as more than a dozen less luminous knots. In the deconvolved 
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near-infrared images, these knots are readily visible. The northeastern knot seems to dominate at 
K' . No additional knots are seen in the near-infrared. 
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